Chemical compass for avian magnetoreception as a quantum coherent device 
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It is known that more than 50 species use the Earth's magnetic field for orientation and naviga- 
tion. Intensive studies particularly behavior experiments with birds provide support for a chemical 
compass based on magnetically sensitive free radical reactions as a source of this sense. However, the 
fundamental question of whether and how quantum coherence plays an essential role in such a chem- 
ical compass model of avian magnetoreception yet remains controversial. Here, we show that the 
essence of the chemical compass model can be understood in analogy to a quantum interferometer 
exploiting quantum coherence. Within the framework of quantum metrology, we quantify quantum 
coherence and demonstrate that it is a resource for chemical magnetoreception. Our results allow 
us to understand and predict how various factors can affect the performance of a chemical compass 
from the unique perspective of quantum coherence assisted metrology. This represents a crucial 
step to affirm avian magnetoreception as an example of quantum biology, which exhibits a direct 
connection between quantum coherence and biological function. 



Despite the growing interest from chemists and biolo- 
gist and more recently researchers from quantum physics 
and quantum information, there remain to date only a 
handful of biological phenomena that are suspected or 
proven to rely on quantum effects. These include impor- 
tant biological processes such as light harvesting [1-6], 
the human sense of smell [7-9] and avian magnetorecep- 
tion [10-15, 17-25]. The effects of weak magnetic field in 
nature have been observed for a long time [26], ranging 
from the growth of plants to the remarkable orientation 
and navigation abilities of birds, mammals, reptiles, am- 
phibians, fish, crustaceans and insects. There are two 
principal mechanisms that may account for avian mag- 
netoreception: one based on magnetically sensitive free 
radical reactions [10-16] and the other on magnetic iron 
minerals [27]. The radical pair mechanism, which was 
put forward to explain a variety of magnetic effects on 
free radical reactions [28] in spin chemistry experiments 
[29], suggests that the avian compass relies on magneti- 
cally sensitive radical pairs formed by photoinduced elec- 
tron transfer reactions [10-15]. The cryptochromes in the 
retina of migratory birds provides a potential physiologi- 
cal implementation of such a mechanism [11, 30-38]. The 
numerous behavior experiments with birds [40], such as 
European robins, have shown that the magnetic orienta- 
tion of birds responds to the inclination of the Earth's 
magnetic field, and is dependent on the wavelength of 
the ambient light [39]. Together with the finding that 
the navigation can be disrupted by very weak external 
oscillating magnetic fields [41], these observations pro- 
vide corroborating evidence for the idea that the chemical 
compass mechanism is involved in avian magnetorecep- 
tion. 

In the last few years, the interest in the long-standing 
problem of avian magnetoreception has quickly extended 
from chemists and biologists to quantum physicists [17- 
25]. Although the detailed biophysical mechanisms un- 



FIG. 1: Quantum interferometer model of chemical 

compass, (a) A quantum system, e.g. a photon or a spin, 
is prepared (BSi) in a coherent superposition of two quan- 
tum states | a) and |/3), which undergoes different evolution 
paths dependent on an unknown physical parameter. The 
measurement after the interference (BS2) of two paths reveals 
the information on the unknown parameter, (b) A chemical 
compass is viewed as a quantum interferometer, the hyperfine 
coupling Hamiltonian Wq recasts the state po of the combined 
system of radical pair and nuclei into a superposition of its 
eigenstates. The weak magnetic field induces the change of co- 
herence phases between different eigenstate evolution paths. 
The spin-dependent chemical reaction channels lead to the 
observable that reveal the interference of coherence phases. 



derlying the chemical compass model of avian magne- 
toreception are quite complex and their details remain 
poorly understood, the radical pair model essentially in- 
volves two unpaired electrons (i.e. a pair of radicals) 
created in spin singlet or triplet states [10, 11, 29]. The 
subsequent dynamics is composed of a quantum coherent 
interaction between the spins as well as an interaction 
with a nuclear spin environment. The former supports 
quantum coherence while the latter is usually considered 
as a noise process that suppresses quantum coherence. In 
previous consideration emphasis was placed on the elec- 
tronic coherence properties and it remained unclear how 
this coherence was related to the function of the chemical 
compass. 

In the present work, we demonstrate that it is essential 



to go beyond this picture and to consider the combined 
system of both radical pair electron spins and the sur- 
rounding nuclear spins as the core of chemical magnetore- 
ception. This leads us to appreciate that the quantum 
dynamics giving rise to the magnetic sensitivity is akin to 
a quantum interferometer in quantum metrology [42-44] . 
Indeed, the global electron-nuclear spin coherence, with 
respect to the eigenbasis of the hyperfine Hamiltonian, is 
found to be essential for the function of a chemical com- 
pass. The validation of this viewpoint depends on an ap- 
propriate quantification of coherence (see supplementary 
information for more discussions). To this end we intro- 
duce a simple and operationally motivated measure of 
quantum coherence, and show its relation with the mag- 
netic sensitivity of a chemical compass. In contrast to 
studies that are only accounting for electron spin coher- 
ence (which we will denote as local quantum coherence 
to differentiate it from the concept of global quantum 
coherence put forward here) [20-24], our approach will 
clearly illustrate the essential role of quantum coherence 
in chemical magnetoreception. 

Quantum interferometer model of chemical compass. — 
In quantum metrology, a system, e.g. a photon or a spin, 
is initially prepared into the following state 

po = |#oX*o|, with l*o) =la\a)+ 7/31/3), (1) 

which is written in the basis of {\a), \/3)} (determined by 
the interferometric element BS1, see Fig.l). A Hamilto- 
nian H(X) = | (|a)(a| — |/3)(/3|) parametrized by a sin- 
gle unknown parameter A results in an evolution that 
changes the relative phase between the two states \a) 
and [/?) to $a = fa — fp, see Fig. 1(a), and leads to the 
state 
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The measurement of the observable M. = |^ r o)(^ r o| ) gives 
the interference (see BS2 in Fig.l) outcome as 

m x (p , * A ) = Tr (p* x M) = 1 - 2| 7q | 2 | 7/3 | 2 (1 - cos$ A ) . 

(3) 
The precision of parameter estimation with such a basic 
quantum interferometer is determined by the contrast of 
the interference fringe as follows 

D 7(Po) = maxm\(p ,$\)-mmm\(po,$\) = 2| 7q | 2 | 73 | 2 . 

(4) 
The role of coherence in such a simple interferometric 
setup can be quantified by 



C = |mA(p§,*A)|* x =o oc |7a| 2 |7/3| S 



with 
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FIG. 2: Effective evolution of a chemical compass un- 
der a weak magnetic field, (a) The change of the abso- 
lute values of the (coherent) density elements (e(t)) averaging 
over time, with e(t) = E m ^„( r mn ~r mn ) 2 ] 1/2 , for the flavin- 
superoxide radical pair with different magnetic field directions 
(6, ip). (b) The singlet yield Ys of the flavin-superoxide rad- 
ical pair as the function of the magnetic field direction (9, ip) 
including (lower) vs. excluding (upper) the coherence phase 

changes <Pmn{t) resulting from the weak magnetic field, see 
Eq.(10). The reaction rates are ks — kr — 0.5/is -1 , and the 
strength of the magnetic field is b — 50/iT. 



which characterizes how coherent the quantum inter- 
ferometer is. The above definition is operational as it 
connects coherence to metrological precision and is thus 
preferable over alternative abstract measures of coher- 
ence, e.g. the summation of the absolute values of the 
off-diagonal density matrix elements (see supplementary 
information for detailed discussions). One can see that 
in the above simple setting of quantum interferometer, 
coherence is an essential ingredient for quantum metrol- 
ogy. Such a framework now allows us to clearly identify 
the role of coherence in chemical magnetoreception ex- 
ploiting the analogy between a quantum interferometer 
and a chemical compass (see supplementary information 
for a complete list of analogies). 

In the standard model of the chemical compass, the 
magnetically sensitive radical pairs formed by photo- 
induced intramolecular electron transfer reactions inter- 
act with a few nearby nuclei via hyperfine couplings, as 
described by the following Hamiltonian 



k=V,A j 



(7) 



where I& . and s^ are the electron and nuclear spin op- 
erators respectively, Tfc denotes the hyperfine coupling 
tensors [29]. For simplicity, we neglect the dipole-dipole 
and exchange interactions between two radicals, which 
is valid when the radical-radical distance is sufficiently 
large or they cancel with each other [45]. The initial 
joint state p of the radical pair and nuclear spins, writ- 
ten in the eigenbasis H^m)} of the hyperfine interaction 
Hamiltonian Hq, can be expressed as 
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In the absence of a magnetic field, the state evolution 



solely under the hyperfine interaction Hamiltonian Ho is 
P n °(t) =J2 r ^ e ' i<n(t) \ u ^^\. (9) 



For a typical radical pair molecule with a few nuclei [46] , 
the leading order effect of the external magnetic field B = 
— gHbb-(sx> + SjCj (with \b\ = 50/iT, which is much smaller 
than the typical interaction strength between the radical 
pair electron spin and the nearby nuclear spins) is to 
introduce an additional magnetic field dependent phase 
on the system coherence, namely 



P 



■H +B(8, V ) 
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E^-^ (t) +^K)KI, 



(to) 

where p w °+ B ( e '^)(i) represents the state evolution un- 
der the total Hamiltonian H = Ho + B(6*, ip) including 
both the hyperfme coupling and the weak magnetic field. 
The effect of the external weak magnetic field might also 
change the absolute values of the coherent off-diagonal 
density matrix elements r mn (m ^ n), which can be 
quantified by e(t) = E m ^ n (^ n - r m „) 2 ] 1/2 - In Fig.2, 
we plot e averaging over time for the possible candidate 
radical pair flavin-superoxide in cryptochrome flavopro- 
tein [14, 15] for different magnetic field directions (0,<p), 
which shows that the change of the absolute values of the 
density matrix elements is very small and has a negligible 
effect on the magnetic sensitivity in the physiologically 
relevant regime (see supplementary information for more 
details). 

The spin-dependent reactions in a chemical compass, 
i.e. the singlet and triplet radical pairs will undergo dif- 
ferent chemical reaction paths and thereby lead to differ- 
ent chemical consequences [11], that witness the magnetic 
field effect on the radical pair dynamics by the magnetic 
anisotropy of reaction yield. For simplicity, we consider 
the scenarios where the reaction rates of the singlet and 
the triplet radicals are identical, i.e. fcg = kx = k (see 
supplementary information for the generalization to more 
general cases). Treating the system dynamics with the 
conventional Haberkorn approach [29], the singlet yield 
can be formulated as 

Y B (Ho,PoA<P)= [ Tr \p n « + ^ e ^(t)M\dt (11) 



with the observable M = ke~ kt (\S){S\ ® kj l dk .). The 
singlet yield in Eq.(ll) can be viewed as a continuous 
generalization of the outcome in a quantum interferom- 
eter as m\(p ,&\), see Eq.(2). With the example of 
flavin-superoxide radical pair, we show in Fig. 2(b) that 
the magnetic anisotropy of singlet yield mainly comes 
from the coherent phase changes ¥?„„(£) induced by the 
external weak magnetic field, as in the interference based 
quantum metrology. 




FIG. 3: The relation between the coherence and the 
magnetic sensitivity of a chemical compass, (a) The 

magnetic field sensitivity Ds("Ho,po) as defined in Eq.(12) 
(upper) and the average value (Ds) (lower) as a function 
of quantum coherence C(rto, po) given by Eq.(14) for 7 x 10 4 
random hyperfme configurations of 5 ~ 6 nuclei with the de- 
polarized initial nuclear spin states. The violet dots repre- 
sent these configurations of reference-and-probe type radical 
pair molecules, which has one radical free of hyperfme cou- 
pling, (b) The magnetic field sensitivity Ds("Hf ,po) as a 
function of coherence measured by C{H^ , po) for the flavin- 
superoxide radical pair (whose hyperfme interaction Hamilto- 
nian is denoted as Ht ) starting from 2 x 10 4 initial states, 
each has the singlet born radical pair and random nuclear 
spin polarizations. In both panels, the reaction rates are 
ks — kr — 0.5/is^ 1 , and the strength of the magnetic field is 
b = 50pT. 



Quantum coherence and magnetic sensitivity. — Fol- 
lowing the inspiration from quantum metrology, we con- 
tinue to quantify the effect of quantum coherence in 
the present scenario of chemical magnetoreception. The 
function of a chemical compass starting from the state po 
can be characterized by the magnetic anisotropy of the 
singlet yield defined as follows 

Ds(%o,po) = ToaaxY s (Ho,po,0,ip) - mmY s (H , Po, 0, <p) 

(12) 
Henceforth we will call Dg('HojPo) as the magnetic sen- 
sitivity of a chemical compass. The coherent part of the 
system state po is denoted as 



QC(H ,po) 






i)("n\ 



(13) 



To quantify the role of coherence in chemical magne- 
toreception, in the same spirit of quantifying the role 
of coherence in the interference based quantum metrol- 
ogy (see Eq.(5)), we introduce the following measure of 
global electron-nuclear quantum coherence for a chemical 
compass with a given hyperfine interaction Hamiltonian 
Ho and a system state po as follows 



C(Ho,Po) = \Ys(n ,gC(H ,Po))\ 



6=0 



(14) 



which represents the contribution of coherence to the sin- 
glet yield in the absence of the external magnetic field. 

To demonstrate the connection between coherence and 
the magnetic sensitivity in a chemical compass, we ran- 
domly sample a large number of hyperfine configura- 
tions and plot the relation between coherence C(Ho,Po) 



and the magnetic sensitivity D^ OjS (yO ), where po — 
\S) (S\ fe . (ld k . /dkj), see Fig.3(a). The result shows that 
typically coherence is necessary for a high magnetic sen- 
sitivity of chemical magnetoreception. Furthermore, we 
also find that the reference-and-probe type of radical 
pairs (namely one radical of which is free from hyper- 
fine coupling) [15, 16] which tends to result in the high- 
est sensitivity usually embody larger coherence. This 
explains such an optimal design principle for a chemi- 
cal compass from the perspective that coherence is a re- 
source for chemical magnetoreception. Considering the 
specific example of the flavin-superoxide radical pair in 
cryptochrome protein, we plot in Fig. 3(b) the magnetic 
sensitivity as a function of quantum coherence from a 
number of randomly chosen initial states of the singlet 
born radical pair but with different nuclear spin polariza- 
tions. The same feature is observed, namely the larger 
coherence the better magnetic sensitivity. We remark 
that the initial state in realistic radical pair reactions 
are usually given by the singlet or triplet radical pair 
state with depolarized nuclear spin states. We assume 
here a larger state space to facilitate the explicit demon- 
stration of how coherence (that changes with the nuclear 
polarization) can directly affect the function of chemical 
magnetoreception. In experiments, this may be achieved 
by chemically induced dynamic nuclear polarization or 
with the assistance of well controllable quantum systems, 
e.g. nitrogen- vacancy centers in diamond [47-49]. In 
fact, both the magnetic field sensitivity and the coherence 
from a initial state with the depolarized nuclear spin state 
are almost identical to the average of these random sam- 
ples with various nuclear polarizations. We also stress 
that the above measure of coherence in Eq.(14) is only 
determined by the hyperfine interaction Hamiltonian Ji 
(without an external magnetic field) and the initial state 
po, yet it well predicts the property (namely the mag- 
netic sensitivity) of the dynamics of a chemical compass 
under a weak external magnetic field. 

An unified picture of decoherence effects. — Following 
the present insight that quantum coherence of the global 
electron-nuclear state is a resource for chemical magne- 
toreception, it is possible to study the effects of differ- 
ent decoherence models on the functioning of a chemical 
compass [21-24] in a unified picture. More specifically, 
we can study how decoherence will destroy the global 
electron-nuclear quantum coherence in a chemical com- 
pass, as characterized by Eq.(14), and thereby deterio- 
rate its magnetic anisotropy of reaction yield. Since a 
chemical compass shall work under ambient conditions, 
the noise from the environment of the core system (i.e. 
the radical pair and the surrounding nuclear spins) will 
inevitably affect its function. The noise effects varies for 
different decoherence models, which can be described by 
the following Lindblad type quantum master equation 
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FIG. 4: The noise effect of difference types of noise on 
the magnetic sensitivity of a chemical compass, (a) 

The magnetic sensitivity Ds(Ho,po) of the flavin-superoxide 
radical pair as a function of the decoherence rate £. (b) The 
magnetic sensitivity Ds(%o, Po) as a function of the coherence 
C(Ho,po) for the local dephasing (red, square), the electron 
spin relaxation (blue, circle), and the singlet-triplet dephasing 
(violet, diamond). The reaction rates are ks — kr — 0.5/us -1 , 
and the strength of the magnetic field is b — 50/j.T. 
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where "H = Ho + B(8, ip) is the total Hamiltonian includ- 
ing the hyperfine interactions Ho and the external mag- 
netic field B(6, (p), Qs and Qt are the projectors into the 
singlet and triplet subspace of the radical pair state indi- 
vidually, and [a;, y]+ = xy + yx. The above master equa- 
tion is based on the Haberkorn approach [29] by adding 
the dissipator C{p) which represent noises. We remark 
that it is not trace preserving as it is restricted in the 
subspace of active radical pair state, and the singlet yield 
can then be calculated as Y s = k s J(S\p(t)\S)dt, which is 
equivalent to Eq.(ll) in the absence of decoherence (see 
supplementary information and Ref. [18]). We consider 
three typical classes of noises, namely the local dephas- 
ing model £/ = {erf ,cr^} with a z = |a)(a| - \0){P\, the 
spin relaxation Cjj = {er±,<7±} with a + = \a)(/3\ and 
(7_ = |/3) (ck|, and the singlet-triplet dephasing model [50] 
Cm = {cgy } with o~st = 2|«S)(«S| — I, where the su- 
perscript D and A represent the two radicals. These 
decoherence models are all relevant in the radical pair 
mechanism but behave very different from each other in 
the sense that they destroy varying aspects of local elec- 
tron spin coherence of radical pairs: they will respectively 
lead to the decay of local electron spin coherence, of elec- 
tron spin longitudinal component, and of singlet-triplet 
coherence at a rate given by £& = £. 

In Fig. 4(a), one can see that the effects of these 
noise models on the magnetic sensitivity of the flavin- 
superoxide radical pair are quantitatively very different 
as a function of the same decoherence rate £. Such a fact 
is natural but gives no general insight into the problem 
how and why different types of noise would affect chemi- 
cal magnetoreception to different extent. Instead, we use 
the proposed measure of coherence C(Ho, Po) and find a 



universal relation between the decoherence effects on the 
magnetic sensitivity and the coherence, which holds for 
different noise models studied here, see Fig. 4(b). This 
supports the observation that the global electron-nuclear 
coherence is indeed responsible for the magnetic sensitiv- 
ity of a chemical compass. Furthermore, by studying the 
effect of noise on the coherence one can predict its effect 
on the magnetic field sensitivity on the basis of its direct 
connection to the global electron-nuclear coherence. 

Discussions. — We have introduced the viewpoint of 
chemical magnetoreception as a quantum interferometer. 
This perspective allows us to reveal that quantum coher- 
ence in the joint basis of the hyperfine interaction Hamil- 
tonian of the radical pair electron spins and the nearby 
nuclear spins is essential for the function of a chemical 
compass. We introduce an operationally defined mea- 
sure for this kind of global electron-nuclear coherence 
to quantify its role in chemical magnetoreception, and 
demonstrate its direct connection with the magnetic sen- 
sitivity of a chemical compass. Our result offers a new 
perspective to understand how a variety of design prin- 
ciples will affect chemical magnetoreception, such as the 
hyperfine configurations, the nuclear spin polarizations 
and the different types of noise. It thus evidences that 
the chemical compass model of avian magnetoreception 
is similar to a coherence based quantum device, and pro- 
vides us with a route to better understand the role of 
quantum effects in chemical magnetoreception. The ver- 
ification of the present observation is conceivable either 
in spin chemistry experiments or by quantum simulation 
with well controllable systems [48] to emulate the radical 
pair dynamics. 
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Supplementary Information 

The hyperfine couplings for the nuclei in the flavin radical. — The flavin-superoxide radical pair is the postulated 
candidate system responsible for avian magnetoreception, the structure of which [SI] is shown in Fig.l. The isotropic 




FIG. 1: The structure of the flavin radical FADH*. The most significant hyperfine interactions in FADH* are from the 
nitrogen N5 and N10, and the proton H5 (blue). 

hyperfine coupling constants (a; so ), the principal values of the anisotropic part of the hyperfine tensors (Ta, with 
i = 1,2,3) and the principal hyperfine axes for nuclei in the avin (FADH*) radical that we use to calculate the 
results for the flavin-superoxide radical pair in both the main text and the following section of this supplementary 
information is based on the data in [SI], see Table I, which include the five nuclei with the largest isotropic hyperne 
couplings and/or the largest anisotropy. 



TABLE I: hyperfine couplings for the nuclei in the flavin radical FADH* 



Nucleus 


a iso (mT) 


Tii(mT) Principal hyperf axes 


N5 


0.393 


0.498 


0.4380 


0.8655 0.2432 






0.492 


0.8981 


0.4097 0.1595 






0.989 


0.0384 


0.2883 0.9568 


H5 


0.769 


0.616 


0.9819 


0.1883 0.0203 






0.168 


0.0348 


0.2850 0.9579 






0.784 


0.1861 


0.9398 0.2864 


N10 


0.212 


0.242 


0.9703 


0.2207 0.0992 






0.234 


0.2383 


0.9426 0.2340 






0.476 


0.0419 


0.2506 0.9672 


Hl'r 


0.390 


0.062 


0.1902 


0.3965 0.8981 






0.033 


0.9156 


0.4017 0.0165 






0.095 


0.3542 


0.8255 0.4395 


H6 


0.158 


0.060 


0.0362 


0.2937 0.9552 






0.044 


0.7948 


0.5879 0.1507 






0.104 


0.6059 


0.7537 0.2546 



The basics of quantum metrology. — In a quantum interferometer, see Fig. la in the main text, the system is prepared 
in a general state written as 



A) = |*o)(*o|, with |^ ) = 7ct |a)+ 7/3 | 



(1) 



when e.g. a photon going through a beam splitter or a f pulse is applied on a spin. The measurement-after-interference 
is given by the observable M. = | v E'o)(*o|, and the interference fringe is 



m\(po, *a) = Tr (p^M) = 1 - 2| 7q | 2 | 7/j | 2 (1 - cos$ A ) 



(2) 



The precision of parameter estimation with such a basic quantum interferometer is determined by the contrast of the 
interference fringe as defined in the following 



D^(po) = maxTO A (p ,$A) - minm A (p , *a) = 2|7 Q | \^p\ = 7. 



(3) 



The contrast of the interference fringe D^(po) illustrates how coherent the quantum interferometer is, as seen in 





FIG. 2 

2|7c| 2 



(a): The interference fringe of a quantum interferometer raj, as a function of the phase $6 with the value of 7 
7,3J 2 = 1 (red, square) and 0.5 (blue, circle), (b): The contrast of the interference fringe D 7 as a function of 7. 



Fig.2(a-b), the more coherence, the larger the interference contrast and thus the better the precision of parameter 
estimation for quantum metrology [S2, S3]. We remark that the quantification of the role of coherence in a quantum 
intcrferometcr-bascd metrology as given in Eq.(4) of the main text is is operational and connects coherence to metro- 
logical precision, and thus is more preferable than the simple summation of the absolute values of the off-diagonal 
density matrix elements, although these two measures are identical for the present observable M. = |^ r o)(^ r o|- 
However, they may not be equivalent for general observables. The key point being that our measure of coherence 
does not merely consider a state but is operationally defined as it takes account of the dynamics of the system (from 
the hyperfine interactions of the radical pair) and the measurement procedure too. Having introduced the basics of 
quantum metrology, we list in Table II the analogies between a quantum interferometer and a chemical compass that 
are relevant for the function of a quantum coherence device, which provide us inspirations to characterize the role of 
coherence in chemical magnetoreception. 



TABLE II: Analogy between a quantum interferometer and chemical compass model 





Quantum Interferometer 


Chemical Compass model 


Basis 


Interference paths 


Eigcnstates of hyperfine interaction Hamiltonian 


Observable 


M = |* >(*o| 


Spin-dependent reactions M — ke~ kt (\S)(S\ ®fc Id k ) 


Outcome 


m A (p ,<I> A ): Eq.(3) 


Singlet yield Y s (Ho,Po, 9, p): Eq.(ll) 


Sensitivity 


Interference fringe contrast D^(po): Eq.(4) 


Magnetic anisotropy of singlet yield T>s{Ho, po)'- Eq.(12) 


Coherence 


C = |ra A (pg,$ A )| $A=0 : Eq.(5) 


C{Uo,po) = \Ys(rl ,gC(rl , P0 ))\ b ^ : Eq.(14) 



The effective evolution of a chemical compass. — The system state of the radical pair and nuclear spins po, written 
in the eigenbasis {|^m)} of the Hamiltonian "Ho, can be expressed as 



Po 



= ^ Tmn \ U m)Wm\, With 



r mn = \v m \p \v n 



(4) 



To show that for a typical radical pair molecule with a few nuclei, where the external magnetic field (much smaller 
than the hyperfine interactions) can be treated as a perturbation, its effect is mainly to introduce an additional phase 
on the system coherence, namely 



p(t) = p^+^^it) « ^r^^e-'KnW+^Wll^X^I, 



(5) 
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FIG. 3: (a) The fidelity (averaging over time) as defined in Eq.(7) between the system state p(t) 



n Ho+B{e, v ) 



(t) and the 



corresponding state neglecting the change of the absolute values of the density matrix elements for the flavin-superoxide radical 
pair with different magnetic field directions (6,f). (b) The fidelity f(t) as a function of time t for the specific magnetic field 
direction (6, </?) = (0, 0). The strength of the magnetic field is b = 50 pT. 



see Eq.(8-10) in the main text, we use the following quantity 



e(t) 




(6) 



to demonstrate that the weak magnetic field would only slightly change the absolute values of the density matrix 
elements when written in the eigenbasis of the hyperfine interaction Hamiltonian "Ho, see Fig. 2(a) in the main text. 
The calculation is done for the flavin-superoxide radical pair in cryptochrome protein, and the hyperfine coupling 
consists of 2 nitrogen (spin-1) and 3 hydrogen (spin-|) nuclear spins, see Table I [SI]. Here, we further support the 
observation by using another quantity to show the slight change in the absolute values of the density matrix elements 
as 



f(p(t),pc(t)) = Tlj y/rif) P (t)yfctf), 



(7) 



which defines the fidelity between the system state p(t) = p H( > +B ( -v) {t^ (under the evolution governed by the total 
Hamiltonian % = Ho+B(6, ip)) and the corresponding state neglecting the change in the absolute values of the density 
matrix elements as 



Pc(t) = 5Z r " 



-*[^»W+A.(*)]| 0<I/fi |. 



(8) 



In Fig. 3, we plot the average value of (f(p(t),pc(t))) for the magnetic field with different directions (9,<p). It can be 
seen that the fidelity is very close to 1, which indicates that the change in the absolute values of the density matrix 
elements is very small. 



Nuclear spin polarization and the magnetic sensitivity of a chemical compass. — In this section we would like 
to compare the magnetic sensitivity with a polarized and unpolarized nuclear environment to provide further 
evidence that the global electron-nuclear coherence is a meaningful coherence measure. To this end, we use 
the flavin-superoxide radical pair as an example and compare the magnetic sensitivity from the depolarized and 
polarized nuclear spin initial states, see Fig. 4. The improvement of the magnetic sensitivity can be attributed to the 
enhancement of the global electron-nuclear spin coherence as measured by C(%o,Po) = |Ys(^o; GCfflo, po))\ b=0 , see 
Eq.(14) in the main text, which increases from 0.105 to 0.2. More generally, the concept of global electron-nuclear 
spin coherence make it possible to understand the effect of nuclear polarization on the magnetic sensitivity of a 
chemical compass straightforwardly. It should be noted that with the techniques of chemical-induced dynamical 
nuclear spin polarization or appropriate quantum control technologies [S4, S5] (partial) polarization of the nuclear 
spins is possible and thus this comparison is also conceivable with current experimental technology. 
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FIG. 4: Nuclear spin polarization and magnetic sensitivity of the fiavin-superoxide radical pair. The dependence of 
the singlet yield Ys(Ho, po, 0, tp) on the magnetic field direction (8, ip) of the fiavin-superoxide radical pair with the depolarized 
(a) and polarized (b) nuclear spin initial state. The common parameters for both panels are: the strength of the magnetic 
field is b — 50/iT, and the reaction rates are ks = kr — 0.5/is" 1 . 



Dynamic features of the the global electron-nuclear spin coherence. — The system state of the radical pair and 
nuclear spins po, written in the eigenbasis {|^ m )} of the Hamiltonian Hq, can be expressed as 



Po 



m.n 



r m n\v m ){vm\, with r. n 



(VmlPoWn) 



In the absence of a magnetic field, the state evolution is as 



(9) 



(10) 



We stress that the definition of the global electron-nuclear spin coherence is only based on the system dynamics 
without a magnetic field, it is thus not a representation of the system dynamics under the action of the weak magnetic 
field (e.g. the Earth's magnetic field). In Fig. 4 of the main text, we demonstrate that the effects of different noises 
can be understood in a unified picture by looking at how each type of noise destroys the global electron-nuclear spin 
coherence. Here we illustrate the decoherence effects in more details by looking at the dynamical features. With a 
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FIG. 5: Upper panels: The dependence of the singlet yield Ys(T-io,po,&, , p) on the magnetic field direction (0,cp) of the 
fiavin-superoxide radical pair under different types of decoherence: without noise (a-1), the local dephasing model (b-1), 
the spin relaxation (c-1), and the singlet-triplet dephasing model (d-1). Lower panels: The singlet fraction (in the state 
evolution p H ° (t) that is only governed by the hyperfine interaction, namely without the magnetic field) as a function of time 
t for different types of decoherence: without noise (a-2), the local dephasing model (b-2), the spin relaxation (c-2), and the 
singlet-triplet dephasing model (d-2). The insets show the zoom in of the time regime [1/US, 3ps]in the red box. The strength 
of the magnetic field is b = 50pT, and the reaction rates are ks = kx = 0.5ps _1 . 
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same decoherence rate £, different noise models disrupt the magnetic sensitivity in varying degrees. In Fig.5(a-l,b- 
l,c-l,d-l), we choose £ = 0.4/xs -1 as an example, and plot the dependence of the singlet yield on the magnetic field 
direction (6, ip) for different types of decoherence. The global electron-nuclear spin coherence is manifested by the 
oscillations of the singlet fraction 



Ps (t) = {S\pK°(t)\S) 






-^"«K|s)(<sk 



(ii) 



The decoherence effects on the oscillations of the singlet fraction (in the state evolution p n ° (£) that is only governed 
by the hyperfine interaction) are also plotted in Fig.5(a-2,b-2,c-2,d-2), which evidence how they destroy the global 
electron-nuclear spin coherence, and thereby decrease the magnetic sensitivity when adding a weak magnetic field. 



Random configurations of hyperfine couplings with nuclei. — The hyperfine coupling tensor with one nuclear spin is 
parametrized by the isotropic hyperfine coupling constants (ai SO ), the principal values of the anisotropic part of the 
hyperfine tensors {Ta, with i = 1; 2; 3) and the three orthogonal principal hyperfine axes (nji, fiz, n k z ) with k = 1; 2; 3, 
e.g. see Table I. In our random sampling, the values of ai SO , Ta and (n k ,fiy,n k ) are generated randomly under the 
realistic constraints: |ai SO | < lmT, |Tn| < lmT, IT22I < lmT, and T33 = T\\ + T22. 



The differentiation between global electron-nuclear spin coherence and local electron-spin coherence. — In the main 
text, we investigate the role of coherence with the concept of global electron-nuclear spin quantum coherence, which 
is essentially different from the conventional definition of local electron-spin quantum coherence. In this section of 
supplementary information, we give more explicit definitions for these two kinds of coherence, and further illustrate 
their contributions to chemical magnetoreception in the next section. The concept of quantum coherence arises from 
the basic principle of quantum mechanics, namely a quantum system can be in a superposition of different orthogonal 
states (called as a set of basis states). Quantum coherence of a given state can be uniquely defined only when a set of 
basis states are specified. Expressed in a specific set of basis, the off-diagonal density matrix elements of a quantum 
state represent its coherent properties. In the case where the radical pair is born in the singlet state and the nuclei 
are depolarized, the system (including the radical pair and the nearby nuclei) state is 



Po = \S){S\(g)(l dk ./d kj ) 



(12) 



where d^ is the dimension of the fcjth nuclear spin, 1^ . is the dk j x dfy identity matrix, and \S) = J\ {\ a P) — |/3o;)) 
is the singlet radical pair state. Written in the conventional local spin basis {\a}, \/3}}, the above state has non-zero 




m 



FIG. 6: The illustration of the state of the combined system of the radical pair and the nearby nuclei. The 

density matrix elements |(po)mn| (with po written in Eq.(13)) (a) in the local electron spin basis and (b) the corresponding 
state, see Eq.(14), by vanishing the local electron spin coherence. For better visibility, the absolute values of the density matrix 
elements are normalized by the largest element. 
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off-diagonal density matrix elements (namely local electron-spin quantum coherence) 



Po 



/o o o\ 
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(13) 



The corresponding incoherent state by vanishing its local quantum coherence, namely the off-diagonal density matrix 
elements as in Eq.(13) is (see Fig. 6a) 



/0 o\ 



lCc( Po ) 



1 
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10 
10 


fc-i 




\0 0/ 





/l o 
1 







(14) 
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The definition of local quantum coherence is thus not dependent on the hyperfine coupling properties of individual 
radical pair molecules, based on which it is not likely to enable us to make a general statement on the role of quantum 
coherence in chemical magnetoreception for different radical pair molecules, see e.g. Ref. [S6]. 

Given an example of radical pair with the hyperfine interaction Hamiltonian Hq, the eigenbasis of which are denoted 
as {|^ m )}j the radical pair state can instead be written in the eigenbasis as follows 



A) 



m.n 



'rnn l^ra/ y^n \ ■> 



ith 



{v m \po\Vn)- 



(15) 



where (I^Vn.)} are the eigenstates of the hyperfine interaction Ho- We call the coherence in this representation as global 
electron-nuclear spin coherence, which is abbreviated as quantum coherence without causing confusion in the main 
text. The coherence properties by this definition depends on both the system state and the hyperfine interactions 
in individual molecules. The concept of global electron- nuclear spin coherence is different from local electron-spin 
coherence. In particular, even the local incoherent state ICc(po), see Eq.(14), may exhibit a certain level of global 
quantum coherence, see an example in Fig. 7. In contrast, the global incoherent state corresponding to a system state 
Po as defined in the following 



lCg(H a ,Po) 



\v m ){v n 



with 



= \Vm\P0\Vi, 



(16) 



does not have global quantum coherence. 





0.6 



m 



FIG. 7: The illustration of the state of the combined system of the radical pair and the nearby nuclei. The 

density matrix elements |(/9o)mn||i/ m },|i/„) (a) in the eigenbasis of the hyperfine interaction Hamiltonian Ho (with three nuclear 
spins that have the largest anisotropy of hyperfine couplings (namely N5, H5, N10, see Table I) in the flavin-superoxide radical 
pair) and (b) the corresponding state, see Eq.(14), by vanishing the local electron spin coherence. For better visibility, the 
absolute values of the density matrix elements are normalized by the largest element. 
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The contributions of global electron-nuclear spin coherence and local electron- spin coherence to the magnetic 
sensitivity. — In the previous studies on the role of quantum coherence (or an extended concept of quantum en- 
tanglement) in a chemical compass, work usually focused on local radical pair electron-spin coherence. In the main 
text, we introduce the notion of global electron- nuclear spin coherence (see also the previous section), and show that 
it is essential for the function of a chemical compass. Here, we compare the contributions of both local and global 
electron-nuclear spin coherence to the magnetic sensitivity of a chemical compass, to provide further support for our 
observation that global electron-nuclear spin coherence is a more fruitful concept allowing us to characterize the role 
of coherence in chemical magnetoreccption. We numerically sample a number of hyperfine coupling configurations 
with from three up to six nuclei, and we assume that the radical pair is born in the singlet so that the initial system 
state is po = |«S)(«S| ®kj (Ifc/rffc)- The contribution of the global electron-nuclear spin coherence to the magnetic 
sensitivity is defined as 



/^ C (D S ) = (D s (Ho,gC(H ,p ))/n s (Ho,p )) = 1 - <D s (Ho,IC e (H ,Po))/D s (H ,Po)), 



(17) 



Where GC(H ,p ) = J2 m7 t n r mnWm}Wn\, lCg(H ,Po) = J2m r mmWm}(^, 

contrast, the contribution of the local electron-spin coherence is 



hi ii with r mn 



(v m \PoWn), see Eq.(16). In 



/"-(D s ) = 1 - (D s (Ho,IC £ (p ))/D s (Ho,Po)>, 



(18) 



where ICc(Po) 1S the local incoherent state, see Eq.(14). 
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FIG. 8: (a) The histogram of the magnetic sensitivity Dg(p) (upper) that is larger than 0.04, and the magnetic sensitivity 
from the corresponding system state with vanishing global electron-nuclear spin coherence Ds(Wo, IC,c(po)) (middle) and with 
vanishing local electron-spin quantum coherence T>s('Ho,ICg(po)) (lower) respectively, (b) The percentage of the magnetic 
sensitivity contributed by the local electron spin coherence and the global electron-nuclear spin coherence quantified by / (Ds) 
and / (Ds) respectively (see the definitions in Eq.(17,18)), averaging over the samples with Ds(po) £ [Ds, Ds + 0.01). In both 
panels, the hyperfine couplings are randomly sampled with in total 3-6 nuclear spins; the reaction rates are k$ — kx = 0.5/xs -1 
and the strength of the magnetic field is b — 50/j.T. 

In Fig. 8(a), we compare the magnetic sensitivity of a chemical compass with the ones excluding local electron-spin 
coherence (namely with the system state JCc(po), see Eq.(14)) and excluding the global electron-nuclear spin 
coherence (namely with the system state ICg(po), see Eq.(16)). It can be seen that the magnetic sensitivity can still 
be very high if the local electron-spin coherence is absent. Instead, the magnetic sensitivity significantly reduces 
when the global electron-nuclear spin coherence is absent. In Fig. 8(b), we plot the contribution of coherence to the 
magnetic sensitivity, as quantified respectively by /^ (Ds), see Eq.(17), and f cc (Ds), see Eq.(18), by averaging over 
the samples with Ds(po) € [Ds, Ds + A). The results show that / ec (Ds) is above 85% and much larger than / £C (Ds), 
which clearly evidences the essential role of the global electron-nuclear spin coherence in chemical magnetoreception. 

Haberkorn master equation for the radical pair dynamics. — To derive the Haberkorn master equation, we introduce 
the notation of |e), \g) to denote the radical pairs when they are excited or in the ground state respectively. In 
the conventional picture, we usually only deal with the excited radical pair state. One can write the system total 
Hamiltonian as follows 



n, 



m 



|e)(e|= H + B(9,cp) 60|e)(e|, 



(19) 
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where Ho represents the hyperfine interactions. We note that the hyperfine coupling and the interaction with the 
magnetic field do not take effect if the molecule is in the ground state \g). The reaction channels from the singlet and 
triplet radical pair are represented by the following Lindblad operators 

Ls = Qs(g)l.9>(e|, (20) 

L T = Qt® ls)(e|, (21) 

where Qs and Qt are the projectors into the singlet and triplet subspace respectively. Therefore the system dynamics 
can described the following master equation as 

d 



,,Peg = -i[Heg,Peg] + ks ( L S p eg L\ - -L^LsPeg - ~Pegp\pS ] + &T ( L T p eg L T - -L T L T p eg - -p eg L T L T 



(22) 



We define the projector into the excite state subspace as 

P e = I0|e)(e| (23) 

and apply it to the master equation in Eq. (22), which leads to 



_d 
dt 

ks 



(P e p eg P e ) = -i[H eg , {PePegPe)} (24) 



., XQS <8> |e)(e|) • (PePegPe) + {PePegPe) ■ (Q S ® |e)(e|)] 

- -f [(Qt (8 |e)(e|) • (P ePeg P e ) + {P e p e gPe) ■ {Qt 2) |e)(e|)] 
We denote P e p eg P e = p e ® |e)(e|, and obtain the Haberkorn master equation for p e from Eq.(25) as follows 

— p e = -i[% Pe] - Y (QsPe + PeQs) ~ "y (QtPc + PcQt) (25) 

where H = Hq + B(#, tp). We remark that the density matrix p e is not trace preserving as it only accounts for the 
population of excited state. The singlet yield (i.e. the singlet ground state population) after time t can be calculated 
as 

Y s (t) = Tr [(Q s ®\g)(g\)p eg (t)}. (26) 

From Eq. (22-25), it is easy to obtain the following dynamics equations for the singlet yield as 

^Y s (t) = k s Tr[Q sPe (t)]. (27) 

that is 

Y s = fc s J Tr [Q sPe (t)]dt. (28) 

In the particular case of ks = &t = k, the solution of p e (t) can be simplified by introducing the density matrix 
pHo+B(e,tp)^^ gee gq .(10) in the main text and Eq.(5), which is the state evolution governed by the system Hamiltonian 
H = H + B{9,<p) as 

Ap*o+S(fl,*0( t ) = -i[H,p n °+^ e <*Xt)}. (29) 

ClTi 

It can be seen from the master equation Eq.(25) that p e (t) — e~ kt p Ho+B ^ ' v '(t). Therefore, the singlet yield in Eq.(28) 
can be obtained as follows (see also Eq.(ll) in the main text) 

Y s (H 0)P(h e,cp)= JTv\p H «+^ ^{t)M\dt, (30) 



with the observable 



M = ke- kt (\S)(S\® ki I dk .), 
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(31) 



which is the so-called exponential model [S7] 



Generalizations to more general radical pair reactions. — In the main text, for simplicity, we have considered the 
scenarios that the reaction rate k$ — &t = k. In such a case, the coherence generated during the evolution is very 
small, see the above section of The effective evolution of a chemical compass. Our idea can also be generalized to 
the radical pair reactions with kg ^ kx by mapping to a continuous type of quantum interferometer so that we can 
also incorporate the coherence generated during the evolution. To this end, we divide the total time evolution t into 
small intervals [At, At, ■ ■ ■ , At], the dynamics of a chemical compass model can be view as a sequence of quantum 
interference processes [ — free evolution for time At — spin dependent chemical reactions — ] N because the singlet 
fraction is continuously leading to reaction products, see Fig. 9. 



B(6,<j>,At) 




FIG. 9: A continuous quantum interferometer model of chemical compass. The dynamics of a general chemical 
compass is viewed as a sequence of quantum interference processes, the hyperfine coupling Hamiltonian "Ho recasts the system 
state p into a superposition of its eigenstates. The weak magnetic field induces the change of coherence phases between different 
paths. The spin-dependent chemical reactions lead to the phase dependent observables. We refer the list of analogies between 
a quantum interferometer and a chemical compass to Table II . 



The evolution of the system is described by the Haberkorn master equation [S7] as 



dt' 



-i[H f p] - -^-[Qs,p}+- —[Q t ,p]+, 



(32) 



where % = 'Hq + ^{9, <p) is the total Hamiltonian including the hyperfine couplings and the external magnetic field, Q$ 
and Qt are the projectors into the singlet and triplet subspace individually, and [a;, y]+ = xy + yx, see also Eq.(25), 
here we simplify the notation p e as p without causing confusion. The measure of coherence, in Eq.(13) of the main 
text, can then be generalized as 



C(H ,Po) = I f k s TY[gC(H , P (t))-M}dt\, with GC{U a ,p{t)) = J^ r mn {t)\v m ){v n \, 



(33) 



with r mn (t) — (v m \p(t)\v n ) . The above quantity characterizes how coherent the corresponding quantum interference 
processes are. The corresponding measure of local electron-spin coherence is defined as 



Cc(Po) 



k s Tr[£C(p(t))-M]dt\, 



(34) 



with CC(p(t)) = p(t) — IC c(p(t)), see the definition in Eq.(14). Using this generalized measure to demonstrate the 
role of coherence in a general chemical compass, we first consider the flavin-superoxide radical pair as an example, 
and vary the ratio kg/k? while fixing the value of fc^. For the simplicity of calculations, we include three nuclear 
that have the largest anisotropy of hyperfine couplings (namely N5, H5, N10, see Table I). In Fig. 10 (ai), it can be 
seen that the magnetic sensitivity Dg changes with the ratio ks/kr, and in the extreme case of ks = 0, obviously 
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D s = 0. The direct relation between the magnetic sensitivity Ds and the global electron-nuclear coherence C, as 
defined in Eq.(33), is shown in the inset of Fig. 10 (ai). For comparison, we also plot in Fig. 10 (ai) the local 
electron coherence Cc, which is negligible even though the magnetic sensitivity can be very high. Similarly, we plot 
in Fig. 10 (02) the magnetic sensitivity excluding the contributions to the singlet yield production (and thereby to 
the magnetic sensitivity) from the local electron spin coherence and the global electron-nuclear spin coherence (as 



sGC 



denoted by D s L and D s 6 respectively). The results show that the global electron- nuclear spin coherence is indeed 
indispensable to the magnetic sensitivity for the present example of flavin-superoxide radical pair. In Fig. 10 (&i,&2), 
we compare the magnetic sensitivity with/without the contributions from the local electron spin coherence and the 
global electron-nuclear spin coherence for a number of random hyperfine configurations. The contributions to the 
magnetic sensitivity by the local electron spin coherence and the global electron-nuclear spin coherence are quantified 
by / £C (D S ) = 1 - D^ Cc /D s and / ec (D s ) = 1 - Dg Cs /D s respectively (similar to the definitions in Eq.(17,18)). The 
results show that / ec (Ds) is above 85% (and even above 90% for the high magnetic sensitivity), which clearly evidences 
the essential role of the global electron- nuclear spin in chemical magnetorcccption. In contrast, the contribution of 
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FIG. 10: (ai): The magnetic sensitivity Ds, the generalized global coherence measure C (Eq.33) and the corresponding measure 
of local coherence Cc (Eq.34) as a function of the ratio ks/kr for the flavin-superoxide radical pair. Inset: The magnetic 
sensitivity Ds as a function of the generalized global coherence measure C as defined in Eq.33). (a2): The magnetic sensitivity 
excluding the contributions from the local electron spin coherence and the global electron-nuclear spin coherence, which are 
quantified by D s c and D s ° respectively, as a function of the ratio ks/kr for the flavin-superoxide radical pair. The common 
parameters for the panels of (ai) and (a 2 ) are: the hyperfine couplings with two nitrogen nuclear spin (with spin 1) and one 
hydrogen nuclear spin (with spin |), namely N5, H5, N10, see Table I. (bi): The histogram of the magnetic sensitivity Ds 
(upper) that is larger than 0.02, and the histograms for the corresponding magnetic sensitivity D s ° (middle) and D s c 
(lower) that exclude the contributions from the global electron-nuclear spin coherence and the local electron spin coherence 
respectively. (b 2 ): The percentage of the magnetic sensitivity contributed by the local electron spin coherence and the global 
electron-nuclear spin coherence quantified by / (Ds) = 1 — D s £ /Ds and / (Ds) = 1 — D s S /Ds respectively (similar to 
the definitions in Eq.(17,18)), averaging over the samples with Ds(po) £ [Ds,Ds + 0.02). In the panels of (bi) and (b 2 ), the 
hyperfine couplings are randomly sampled with one radical coupled with three nuclei and the other radical is free of hyperfine 
interaction, the ratio between the singlet and the triplet reaction rate ks/kr is randomly chosen in the regime [0, 1]. The total 
number of random samples is 1.2 x 10 4 . The common parameters for all the four panels are: the strength of the magnetic field 
is b — 50fiT, and the triplet reaction rate is kr = 1/is -1 , the initial system state is po = \S) {S\ ® fe ( Id k /dk ) ■ 



the local electron spin coherence (see / (Ds) in Fig. 10 (62)) is much less significant. 
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FIG. 11: (ai): The magnetic sensitivity Ds, the generalized global coherence measure C (Eq.33) and the corresponding local 
coherence Cc (Eq.34) as a function of the ratio ks/kx for the flavin-superoxide radical pair, (aa): The magnetic sensitivity 
excluding the contributions from the local electron spin coherence and the global electron-nuclear spin coherence, which are 
quantified by D g c and D s 9 respectively, as a function of the ratio ks/kx for the flavin-superoxide radical pair. In the panels 
of (ai) and (a2): the hyperfine couplings consist of two nitrogen nuclear spin (with spin 1) and one hydrogen nuclear spin (with 
spin |), namely N5, H5, N10, see Table I. (bi): The comparison between the magnetic sensitivity Ds and the corresponding 
one D s g excluding the contributions from the global electron- nuclear spin coherence. (b2): The magnetic sensitivity Ds as 
a function of the generalized global coherence measure C (Eq.33) as defined in Eq.(33). In the panels of (bi) and (b 2 ), the 
hyperfine couplings are randomly sampled with one radical coupled with three nuclei and the other radical is free of hyperfine 
coupling, the ratio between the singlet and the triplet reaction rate ks/kT is randomly chosen in the regime [0, 1]. The total 
number of random samples is 2 x 10 4 . The common parameters for all the four panels are: the strength of the magnetic 
field is b — 50/uT, and the triplet reaction rate is kx — Ifis' 1 , the initial system state is chosen as the hypothetical state 

Po = (l4/4)(g) fe . (ld kj /d hi ) [S8]. 



In Ref. [S8], it was shown that in the case of ks ^ kx, a chemical compass starting from the hypothetical initial 
state pa — QLi/4.) (££) fe (l^./dkA, which has no coherence expressed in any set of basis, can still have a certain 
level of magnetic anisotropy in the singlet yield. Although such an initial system state is quite unusual in realistic 
radical pair reactions, it is an interesting hypothetical scenario to discuss about the role of coherence in chemical 
magnetoreception. It is crucial to note that the imbalance of singlet and triplet reaction channels will introduce 
global coherence [S8], which can be captured by the proposed measure in Eq.(33) as we will show in the following. 
In Fig.ll(ai,a2), we take the flavin-superoxide radical pair as an example, and vary the ratio kg/kx while fixing 
the triplet reaction rate kx- The results demonstrate the similar essential role of the global electron- nuclear spin 
coherence in this scenario. In Fig. 11 (bi, 62), we calculate the magnetic sensitivity Dg for a number of radical pair 
molecules with random hyperfine configurations, in which one radical is coupled with three nuclei while the other 
radical pair is free of hyperfine interaction, and the ratio kg/kx is randomly distributed in [0, 1]. It can be seen that 
the magnetic sensitivity D s is likely to be very small, nevertheless, one can still witness a direct connection between 
the global electron-nuclear spin coherence and the magnetic field sensitivity. 
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